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ABSTRACT 

We present nonlinear theoretical results concerning the predicted pulsational behaviour of 
stellar models suitable for the sample of Cepheids belonging to the Large Magellanic Cloud 
(LMC) star cluster NGC 1866. The blue and red edges of the instability strip, transformed 
into the observative plane (V, B-V) by adopting current values for the distance modulus 
(DM=18.57 mag., Welch et al. 1991) and for the reddening (E(B-V)=0.06 mag., Arp 1967; 
Walker 1974; Walker 1987), are in agreement with the observed distribution of Cepheids in 
the quoted cluster. Moreover, the distribution of observational data in the Bailey diagram (V 
amplitude versus period) appears in agreement with theoretical predictions derived adopting 
the luminosity level predicted by canonical evolutionary models. The role played by the 
direct ah initio integration of the coupling between pulsation and convection and by the 
input physics on the smoothness of light and velocity curves is briefly discussed. 
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1 INTRODUCTION 

The LMC globular cluster NGC 1866 has been already used 
extensively for testing stellar evolutionary models and for 
shedding light on some problems currently debated in the 
literature (Chiosi et al. 1989; Brocato, Castellani & Pier- 
simoni 1994). This cluster is also expected to play a key 
role in the theory of stellar pulsation, since it provides a 
beautiful sample of 23 Cepheids with a common value for 
both distance modulus and reddening. Welch et al. (1991, 
hereinafter referred to as WCFMM) already disclosed that 
this sample of radial pulsators shows a well defined instabil- 
ity strip which is not revealed by current data on Galactic 
Cepheids, and represents a fundamental benchmark for both 
modal stabilities (fundamental and first overtone) and pul- 
sational amplitudes predicted by theoretical models. 
As a first step of an extended and detailed investigation of 
Cepheid hydrodynamical models we use NGC 1866 to test 
both the accuracy and the internal consistency of the theo- 
retical framework we developed for evaluating the approach 
to limit cycle stability and the predicted pulsational ampli- 
tudes. The main aim of this paper is to discuss the com- 
parison between a new set of limiting amplitude, nonlinear, 
nonlocal and time-dependent convective models and the ob- 
servational constraints provided by this cluster. In §2 the 
theoretical results are briefly discussed while §3 presents the 
comparison with the sample of Cepheids in NGC 1866 for 
which magnitudes, colours, and pulsational amplitudes are 
currently available. The last section summarizes the main 
results of this investigation and outlines future theoretical 
and observational perspectives. 



2 THEORETICAL FRAMEWORK 

According to current evolutionary predictions for Cepheids 
in NGC 1866, we adopted a fixed stellar mass (M/M© =5.0). 
Even though no direct estimation of the metallicity of this 
cluster has been provided yet, we adopted the metal abun- 
dance obtained by Russell & Bessell (1989) and Russel & 
Dopita (1990) for young LMC stars, i.e. Z=0.008. We 
also adopted the helium content Y=0.25 provided by Du- 
four (1984) on the basis of He abundance evaluations in H 
II regions. Similar physical parameters have also been as- 
sumed by Chiosi et al. (1992,1993 and references therein) in 
their extensive investigation of the pulsational properties of 
Cepheids linear models. For properly evaluating the topol- 
ogy of the instability strip we computed four sequences of 
both linear and nonlinear models located at different lumi- 
nosity levels (log L/Lq = 2.80, 3.07, 3.3, 3.6) which cover a 
wide effective temperature range (4800 <Te < 6700 K). The 
location of the blue and red edges of fundamental and first 
overtone pulsators has been evaluated by adopting, close to 
these boundaries, a temperature step of 100 K. The theoret- 
ical framework adopted for evaluating the pulsational prop- 
erties of Cepheids has already been described in a series of 
previous papers (Bono & Stellingwerf 1993,1994; Bono et 
al. 1997b and references therein). In this section we only 
discuss the main differences with the quoted investigations 
concerning the physical assumptions adopted to approach 
Cepheids pulsational models. 

For providing a good spatial resolution throughout the en- 
velope model, the static structures have been computed by 
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adopting an outer boundary optical depth r = 0.0001. At 
the same time, the mass ratio (h) between consecutive zones 
has been assumed equal to 1.04 in the stellar layers located 
at temperatures lower than 6.0 x 10^ K, whereas for higher 
temperatures we adopted h=1.2. The inner boundary has 
been fixed so that the base of the envelope was approxi- 
mately located at a distance from the stellar centre of the 
order of 10% of the photospheric radius (i.e. ro = 0.1 xil^^). 
Unlike envelope models of low-mass stars, these assumptions 
ensure that the mass included in a typical Cepheid model 
ranges from 40 to 50 percent of the total stellar mass. Both 
linear and nonlinear models have been computed by adopt- 
ing the radiative opacity tables recently provided by Iglesias 
& Rogers (1996) for temperatures higher than 10,000 K, 
and the molecular opacities provided by Alexander & Fer- 
guson (1994) for lower temperatures. The method adopted 
for handling the opacity tables is described in Bono, Incerpi 
& Marconi (1996). 

The nonlinear analysis has been performed by perturbing 
the linear radial eigenfunction of the first two modes with a 
constant velocity amplitude of 10 km/sec. In contrast with 
nonlinear RR Lyrae and Type II Cepheid models, which 
have been computed with timesteps fixed by the CFL condi- 
tion for properly handling the development and propagation 
of the shock front along the pulsation cycle, the nonlinear 
Cepheid models have been computed by adopting a number 
of timesteps per period which range from 300 for first over- 
tone models to 400 for fundamental pulsators. This choice 
is due to the evidence that along a full cycle Cepheid models 
develop only mild shocks in the hydrogen and helium ion- 
ization regions. As a test, few selected models computed by 
adopting the CFL conditions do not present any substantial 
difference in the pulsation amplitudes in comparison with 
the models computed by adopting a fixed timestep. The 
number of pulsation cycles needed to approach the asymp- 
totic amplitudes depends on the location of the model in- 
side the instability strip and ranges from 1,000 to more than 
10,000. 

A thorough comparison of our models with both linear 
(Chiosi et al. 1993; Morgan & Welch 1996) and nonlin- 
ear radiative Cepheid models (Sebo & Wood 1995; Simon & 
Kanbur 1995; Buchler et al. 1996), together with a detailed 
analysis of the physical structure of the envelope models will 
be discussed in a forthcoming paper (Bono & Marconi 1997). 
The top panel of Figure 1 shows a selected sample of fun- 
damental and first overtone bolometric light curves located 
at log L/Lq = 3.07, i.e. at the luminosity level predicted 
by canonical evolutionary models. The first interesting re- 
sult is that moving from the blue to the red edge of the 
instability strip the light curves show smooth changes along 
the pulsation cycle and do not present, at least in this pe- 
riod range, any peculiar bump and/or dip. The comparison 
with results available in the literature (Karp 1975; Carson & 
Stothers 1984) strengthens the key role played by the inclu- 
sion of a time-dependent treatment of convective transport 
which reduces the pulsational amplitudes and at the same 
time ensures a smooth excursion of physical variables across 
the ionization regions. The bottom panel of Figure 1 shows 
the surface radial velocities of the same models shown in the 
top panel. The smoothness of these curves over the pulsation 



Figure 1. Top: Selected bolometric light curves as a func- 
tion of the pulsational phase. Solid and dashed lines are 
referred to first overtone and fundamental pulsators re- 
spectively. The magnitude scale is referred to the top 
curve; the other ones have been artificially shifted by 
-0.60 mag. The nonlinear periods are also labelled. Bot- 
tom: Surface radial velocity curves of the models plotted 
in the top panel. The velocity scale is referred to the top 
curve; the other ones have been artificially shifted by -35 
km/sec. The effective temperatures are also labelled. 

cycle is the result of a long ab initio integration time during 
which the radial motions approach their asymptotic ampli- 
tudes and the spurious high-order modes introduced by the 
initial perturbation settle down. A more detailed sum- 
mary of the nonlinear observables obtained in the present 
investigation for both fundamental and first overtone pul- 
sators is given in Table 1. This table presents from left to 
right the logarithmic luminosity, the effective temperature, 
the nonlinear pulsation period, and the fractional radius 
variation, i.e. AR/Rph = (^^"""^ - R"'''")/Rph, where Rp,, 
is the photospheric radius. Column (5) reports the "pure" 
radial velocity amplitude of the surface zone, i.e. no correc- 
tion factor for limb darkening has been taken into account. 
Column (6) gives the bolometric amplitude, while the two 
subsequent columns report the amplitude of static (7) and 
effective (8) logarithmic surface gravity. The static gravity 
is referred to the surface zone, whereas the effective gravity 
has been evaluated by adopting the method suggested by 
Bono, Caputo & Stelhngwerf (1994). Columns (9) and (10) 
give the amplitude of the surface temperature. The former is 
referred to the temperature of the outer boundary, whereas 
the latter is derived taking into account the surface varia- 
tions over a full cycle of both luminosity and radius. The 
final column shows the pulsational amplitude in the V band 
(sec §3 concerning the transformation into the observative 
plane) . 

In order to point out who is doing what concerning the driv- 
ing and the damping of the pulsation instability. Figure 2 

shows the total work curves referred to a fundamental model 
centrally located in the instability strip (log L/Lq = 3.07, 
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Te=5800 K). The light and velocity curves of this model 
arc shown in Figure 1. Since we are interested in disclosing 
the destabilization effects of different elements (H, He, met- 
als) , the total work terms have been evaluated according to 
the relation dW/d{log < T >), where < T > is the time 
average temperature per zone over a full cycle (see Bono 
et al. 1997a). Note that in the literature? the total works 
are generally evaluated as a function of the exterior mass 
(i.e. dW/d{logMe)). According to Figure 2 we find that H 
and He, as expected, are the main destabilization sources 
whereas the Z-bump, at least for the assumed metallicity, 
barely affects the pulsation stability. On the other hand, 
the region located between the HeH ionization zone and the 
Z-bump is the main radiative damping region of the enve- 
lope. Moreover, the artificial viscosity pressure does not 
affect at all the physical structure of this model, since it 
is vanishing throughout the envelope. Turbulent and eddy 
viscosity pressure terms supply a non-negligible amount of 
dissipation in the hydrogen and helium ionization regions 
which in turn reduces the pulsational amplitudes. 

In order to easily identify the opacity bumps which 
cause the appearance of both driving and damping regions, 
the bottom panel of Figure 2 shows the opacity excursions 
throughout the envelope structure. The opacity bump 
which appears in the outermost regions is worth mention- 
ing. This opacity increment, at densities typical of Cepheid 
surface layers, is mainly due to H~ and partially to the 
Rayleigh scattering of hydrogen (see Fig. 7 in Alexander & 
Ferguson 1994). In fact, during a large part of the pulsation 
cycle the surface layers of this model oscillate at tempera- 
tures lower than 5,000 K, reaching the maximum tempera- 
ture excursion (T=5,250 K) only around the phases of the 
luminosity maximum. 

This feature, which has a negligible effect on the pulsation 
driving, could play a key role in the smoothness of both 
light and velocity curves. In fact, thanks to the inclusion of 
molecular opacities, the radial displacement of surface layers 
located above the hydrogen ionization region is governed by 
the physical excursion of temperature and density. 

3 COMPARISON WITH OBSERVATIONS 

In this section theoretical results arc compared with the pho- 
tometric data for NGC 1866 provided by WCFMM. These 
data refer to only nine out of the 23 Cepheids identified in 
NGC 1866, since photometry is not provided in the quoted 
paper for the remainder. In a subsequent paper Welch & 
Stetson (1993) applied their variable star detection tech- 
nique to several frames centered on NGC 1866 and provided 
periods for all but one variable of the whole sample. How- 
ever, both mean magnitudes and colours presented in that 
paper are affected by a poor photometric accuracy due to 
the crowding of central regions and thus are of no use. 
Theoretical boundaries for both fundamental and first over- 
tone instability strips (Z=0.008, Y=0.25 and M = 5Mq) 
have been transformed into the observative plane V-(B- 
V) by using Kurucz's bolometric corrections and colour- 
temperature relations (Kurucz 1992). In Figure 3 the trans- 
formed boundaries are shown together with the observed 
Cepheids of NGC 1866 studied by WCFMM. Following 
WCFMM discussion we have assumed for the cluster a red- 
dening value E(B-V)=0.06 mag. and a distance modulus 



Figure 2. Top: Full amplitude nonlinear total work 
curves versus the logarithmic time average temperature, 
surface at right. This fundamental model is centrally 
located in the instability strip (logL/ Lq=3.07; Te=5800 
K). The solid line shows the total work, the dashed line 
the turbulent work (i.e. the work due to both turbulent 
and eddy viscosity pressure), and the dashed-dotted line 
the artificial viscosity work. Positive and negative areas 
denote driving and damping regions respectively. Bot- 
tom: Logarithmic opacity of the envelope as a function 
of the logarithmic time average temperature. The opac- 
ity profile has been plotted at each time step over a full 
cycle. The arrows mark the location of the main opacity 
sources. 

DM=18.57 mag. Solid lines represent, left to right, the 
First Overtone Blue Edge (FOBE), the Fundamental Blue 
Edge (FBE), the First Overtone Red Edge (FORE) and the 
Fundamental Red Edge (FRE). The FOBE and the FBE 
intersect at log L/Lq « 3.4 and therefore we expect that 
above this luminosity level the instability strip is populated 
only by fundamental pulsators. Indeed at log L/Lq = 3.6 
the nonlinear models computed by adopting the linear first 
overtone radial eigenfunction do not show a stable limit cy- 
cle. On the other hand, toward lower luminosities, the FBE 
and the FORE intersect at log L/L0 = 2.8, thus removing 
completely the "OR region", i.e. the region of the instabil- 
ity strip where pulsators arc characterized by a stable limit 
cycle in the fundamental and in the first overtone mode. 

The agreement with observational data appears satisfactory, 
since most of the observed variables are located well inside 
the theoretical instability strip around a magnitude level 
corresponding roughly to log L/Lq = 3.07. In agreement 
with early results obtained by Bertelli et al. (1993), the 
comparison also suggests the occurrence of first overtone 
pulsators among Cepheids in NGC 1866. The only deviant 
variable is HV 12204, already pointed out by WCFMM for 
its peculiar brightness and blueness. However, from the 
analysis of radial velocities these authors concluded that HV 
12204 is not a cluster member. Two further features of the 
instability strip are worth mentioning: 
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Figure 3. Fundamental and first overtone instability 
boundaries for Z=0.008, Y=0.25 and M = 5Mq. The ob- 
servational data are referred to Cepheids in NGC 1866. 
See text for further details. 

1) moving from higher to lower luminosities the instar- 

bility strip becomes narrower; in the investigated luminosity 
range the colour width roughly decreases of a factor of two 
as the outer edges of the instability region present different 
slopes. These findings support the results obtained by Fer- 
nie (1990) for Galactic Cepheids and the preliminary trend 
suggested by Olszewski (1995) for Cepheids in LMC clusters. 

2) The region of the instability strip in which only first 
overtones show a stable nonlinear limit cycle becomes larger 
at lower luminosities. This outcome supplies a straightfor- 
ward but qualitative explanation for the large amount of first 
overtone pulsators -approximately 30 % of the total sample- 
recently detected by microlcnsing experiments (Boaulieu et 
al. 1995; Cook ct al. 1995) and supports the hypothesis 
originally suggested by Bohm-Vitense (1988) concerning the 
occurrence of such variables among short-period Cepheids. 

In Figure 4 the theoretical relations between periods and V 
amplitudes for both fundamental and first overtone models 
and for the labelled assumptions on luminosity are com- 
pared with the observational data provided by WCFMM. 
Bolomotric amplitudes have been transformed into the V 
Johnson band by adopting Kurucz's static atmosphere mod- 
els. The V magnitude curve can be used to derive either 
a magnitude-weighted amplitude or an intensity-weighted 
one, which is subsequently converted into magnitude units. 
The differences between amplitudes obtained from these two 
different approaches turned out to be negligible (less than 
0.001 mag.). 

It is worth noting that while the comparison in the HR dia- 
gram (instability strips) could be affected by uncertainties in 
the bolometric corrections, in the colour-temperature rela- 
tions, in the reddening, and in the distance modulus, as well 
as in the method adopted to average the observed colours 
and magnitudes over the pulsation cycle, the comparison 
in the Ay-log P plane is independent of almost all these 



Figure 4. Comparison in the Bailey diagram between 
predicted and observed V amplitudes and periods. Sym- 
bols are referred to different luminosity levels (see la- 
belled values). Fundamental and first overtone ampli- 
tudes are plotted by adopting solid and dashed lines re- 
spectively. Asterisks show the observed V amplitudes 
provided by WCFMM. 

troublesome effects, and is therefore far more robust. As 
a consequence, we can infer that Cepheids in NGC 1866, 
at least those for which photometric data are given in the 
literature (WCFMM), are well reproduced in the Bailey di- 
agram by a sequence of pulsational models with Z=0.008, 
Y=0.25, M = 5Mq and log L/Lq ~ 3.07, in remarkable 
agreement with the adopted evolutionary predictions. Note 
that the sequences of models characterized by the same val- 
ues of stellar mass and chemical composition but by luminos- 
ity levels different from log L/Lq = 3.07 are clearly unable 
to match observational data. This means that, in principle, 
the comparison in the period-amplitude diagram proves to 
be a powerful instrument for estimating the luminosity level 
of cluster Cepheids. 

A further interesting aspect worth being investigated in 
order to provide both a suitable test of the present the- 
oretical scenario and a tight constraint on the pulsation 
characteristics of classical Cepheids is the comparison be- 
tween the Fourier parameters of theoretical and observed 
light curves (Simon & Kanbur 1995). Even though such 
analysis is clearly beyond the scope of this paper, a compre- 
hensive comparison with observational data and the impact 
of the pulsation/convection interaction on observables will 
be discussed in a forthcoming paper (Bono & Marconi 1997). 

4 CONCLUSIONS 

The agreement between Cepheid theoretical models and 
the photometric data available in the literature provides a 
sound support to the plausibility of the physical assump- 
tions adopted in the development of this new pulsational 
scenario. In a homogeneous theoretical context wo derived 
both the blue and the red edges of the instability strip. In 
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particular, it is worth noting that the red edges have been 
evaluated without invoking any ad hoc assumption on the 
efficiency of the convective transport over the pulsation cy- 
cle. At the same time present limiting amplitude calcula- 
tions strongly support the existence of first overtone pul- 
sators among Cepheids. In fact, we found that there is a 
well defined region of the instability strip in which only first 
overtone pulsators present a nonlinear limit cycle stability. 
On the other hand, taking into account the distribution of 
pulsational amplitudes in the Bailey diagram we find that 
the luminosity level of Cepheids along the blue loop is in 
agreement with the value predicted by canonical evolution- 
ary models. 

However, a deeper insight on the matter is prevented by 
the limited sample of Cepheids currently available. What is 
relatively surprising about photometric data of young LMC 
clusters (see Table 2 in Welch, Mateo & Olszewski 1993 and 
references therein) with a good Cepheid sample is that even 
though they are the keystone for determining the P-L and 
the P-L-C relations on which rely the distance evaluations 
of the Local Group Galaxies, and in turn the calibration of 
secondary distance indicators, we still lack a comprehensive 
observational scenario for these objects. 
Moreover, although several approaches have been recently 
undertaken for evaluating the metal abundance of stellar 
clusters in Magellanic Clouds (Geisler & Mateo 1989; Ol- 
szewski et al. 1991), a detailed analysis of the chemical 
composition spread of Cepheids in Galactic and Magellanic 
Clouds clusters is even more urgent. 

As far as future theoretical developments are concerned, we 
plan to extend the present nonlinear pulsational scenario 
by taking into account different stellar masses and chemi- 
cal compositions for properly disclosing the dependence of 
Cepheids behaviour on these astrophysical parameters. 
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65.74 


1.078 


0.10 


0.69 


1350 


1700 


1.262 


3.07 


6000 


3.0360 


0.128 


70.40 


0.938 


0.11 


0.74 


1150 


1450 


1.129 


» ■ 


3.07 


5900 


3.2112 


0.126 


68.11 


0.725 


0.11 


0.73 


950 


1200 


0.894 


■ 


3.07 


5800 


3.3979 


0.115 


59.99 


0.624 


0.10 


0.63 


850 


1050 


0.785 


> : 


3.07 


5700 


3.6030 


0.095 


47.62 


0.478 


0.08 


0.49 


650 


800 


0.600 




3.07 


5600 


3.8147 


0.053 


25.43 


0.225 


0.04 


0.25 


300 


400 


0.295 


^; 


3.30 


6000 


4.7509 


0.084 


37.53 


0.712 


0.08 


0.47 


950 


1200 


0.956 




3.30 


5600 


6.0050 


0.086 


38.09 


0.444 


0.08 


0.42 


550 


700 


0.582 




3.30 


5500 


6.3743 


0.068 


29.98 


0.317 


0.06 


0.32 


400 


500 


0.429 




3.30 


5400 


6.7647 


0.043 


17.54 


0.183 


0.04 


0.19 


250 


300 


0.254 




3.60 


5700 


10.3780 


0.129 


45.57 


0.877 


0.11 


0.56 


1100 


1350 


1.123 




3.60 


5600 


11.0760 


0.176 


61.51 


1.186 


0.15 


0.86 


1500 


1850 


1.515 




3.60 


5300 


13.4600 


0.212 


66.52 


1.033 


0.18 


0.90 


1300 


1600 


1.451 




3.60 


5100 


15.1911 


0.146 


40.11 


0.585 


0.13 


0.57 


800 


950 


0.870 




3.60 


5000 


16.2010 


0.046 


11.98 


0.136 


0.04 


0.16 


200 


250 


0.217 



Logarithmic luminosity (solar units). Effective temperature (K). Period (days). ^ Fractional radius variation. 
Radial velocity amplitude (Km/s). ^ Bolometric amplitude (mag.). ^ Amplitude of logarithmic static gravity. 
Amplitude of logarithmic effective gravity. * Surface temperature variation. ^ Effective temperature variation. 
Visual amplitude (mag.). 



